Efficient creation of maximally entangled states by modulation of tunneling rates 
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For systems described by the two-site Bose-Hubbard Hamiltonian, we show that a sinusoidal 
modulation of the tunneling matrix element assists higher-order co-tunneling processes. Using this 
mechanism, we propose an efficient new scheme for creating a coherent superposition of states in 
which all particles are either on one site or all on the other site, the so-called NOON state. This 
scheme yields an almost perfect NOON state periodically. For larger numbers of particles, further 
reduction of the time to create the state is possible if more than one modulation frequency is 
employed. With this scheme, NOON states with a larger number of particles could be realized with 
state-of-the-art techniques for cold Bose gases in a double-well potential. 

PACS numbers: 03.75.Lm, 03.75.Gg, 42.50.Dv, 03.65.Xp 



One of the current challenges in quantum control is to 
create particular sorts of entangled states that are a su- 
perposition of states that are maximally different. The 
largest number of particles for which such states can be 
engineered is six for ions in traps [l[ and ~ 10 photons 
for photons. Atoms in a double- well potential are a 
fundamental system for quantum-state engineering and it 
is a promising candidate for creating such superposition 
states. Especially the coherent superposition of a state 
in which all N particles are in the right well and that in 
which all particles are in the left well (a so-called NOON 
state [H ) is a maximally entangled state 0, Q , and it is 
an important resource for matter- wave interferometry: it 
provides the ultimate quantum limit of the phase reso- 
lution ~ l/Af) i- e -> the Heisenberg limit, for any N Q 
(see also yy). However, NOON states with larger N are 
fragile with respect to decoherence and thus it is crucial 
to create them in a short time. 

In Ref. we studied the tunneling of bosons in a 
double-well potential M d escribed by the two-site Bose- 
Hubbard Hamiltonian [lOj : 



H 



-J(c R c L + c\c R ) + — (c R c R c R c R + c\c L c L c L ), (1) 

where c R and create bosons in the right and left 
well, respectively, J is the tunneling matrix element, 
and U is the on-site interaction. There, we found that, 
starting from a situation in which all particles are in 
the right or the left well, a NOON state is formed by 
co-tunneling provided the interaction is strong enough: 
k/2 = UN/2J ^> 1. In this situation, single-particle tun- 
neling is suppressed by the energy conservation and all 
N particles are forced to be in the same well. While the 
classical theory predicts self-trapping [IJ [lj| , a superpo- 
sition of all particles in the right and left wells is allowed 
in quantum mechanics and an oscillation between these 
two states occurs due to higher-order co-tunneling, which 
yields a NOON state after a quarter of the oscillation pe- 
riod @|. Indeed, a NOON state for N = 2 would have 
been realized by this mechanism in experiments of Ref. 
[HI . A great advantage of this mechanism is that almost 



perfect NOON states are obtained at periodic intervals in 
double- well potentials unlike, e.g ., the well-known proto- 
col of Ref. [l4j (see also [HI, [l6[ and references therein; 
protocols for different setups are described in, e.g., Refs. 
USUI). On the other hand, a serious difficulty with the 
mechanism of Ref. Q is that the time to form NOON 
states is much larger than for other protocols and in- 
creases exponentially with N. In the present work, we 
shall show how this drawback may be overcome. We also 
show that, using our present scheme, NOON states with 
N comparable to the largest value attained for trapped 
ions [l| can be realized with cold atoms under state-of- 
the-art experimental conditions 13]. 

Let us first evaluate the energy splitting AEan of the 
two degenerate states with the same value of |AiV| = 
\Nr — Nl\, where N R and Nl are the numbers of par- 
ticles in the right and left wells, respectively. The on- 
site interaction term, which reads H = U(c R c R c R c R + 
c[c[c L c L )/2 = UAN 2 /4 + U{N 2 - 2A>)/4, gives the 



zcroth-order result for the energy as E y ^' N = UAN 2 /A 



(o) _ 

(see Fig. [TJ. Here N = c R c R + c< L c L and AN = 

c R c R — c l cl. The splitting AEan can be calculated by 
| A7V|th-order perturbation theory treating the hopping 
term as a perturbation and we obtain 



AE 



AN 



2JN K ^ AN ^ 1 f3(N, AN), 



with 



[3(N,AN) = 



[(N+\AN\)/2}\ N 



AJVI-2 



[(JV-|AiV|)/2]! [(|AA|-1)! 



(2) 



(3) 



The tunneling period Tajv between these states is given 
by Tan = 2irfi/ AEan ■ Observe that, for k > 1, Tan 
grows exponentially with |AiV|. 

Our basic idea for reducing the formation time of the 
NOON state is to modulate the tunneling matrix element 
we make a resonance between the states of |AiV| = N 
and | A AT | = TV — 2 (see Fig. [TJ, and take advantage of 
the tunneling between the two states of \AN\ = N — 2 
whose period 2V-2 is much shorter than Tn- Using 
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FIG. 1: (Color online) Schematic diagram of the present pro- 
tocol for reducing the formation time of the NOON state. 



Eqs. @ and ([3]), we estimate the reduction ratio as 
T/v/Tjv— 2 = K 2 [(N - 1)(N - 2)] 2 /JV 3 ~ k 2 N. For, e.g., 
JV = 5 and U/J = 4, the reduction ratio is Tjv/Tjv-2 = 
460.8. From the general point of view, this mechanism 
in which the higher-order co-tunneling is assisted by the 
time-dependent tunneling matrix element would also be 
relevant to other contexts such as in Josephson qubits 
(e.g., [19] for review). Moreover, our theory directly ap- 
plies also to another situation: cold Bose gases in two 
hyperfine states coupled by an external pumping field 
(see, e.g., j7~> ). In this case, J corresponds to the Rabi 
frequency of the external field and by controlling its in- 
tensity one can change J. 

We consider sinusoidal modulation of the tunneling 
matrix clement and employ the two-site Bosc-Hubbard 
Hamiltonian ((I) with J replaced by J(t) where 



J(t) = J(l + Asmcut). 



(4) 



Here A and uj are the amplitude and the frequency of 
the modulation, respectively. This modulation can be 
realized by controlling the height of the barrier in the 
case of the double- well potential (see, e.g., Ref. [2GJ for 
reproducing J(t) of Eq. (HJ) in three-dimensional double- 
well potentials). The condition for a resonance between 
the states with \AN\ = JV and |AJV| = JV - 2 is 



huj~E 



(0) 



E 



(o) _ 

N-2 ~ 



U(N-l). 



(5) 



As we shall see later, to obtain the periodic dynamics, 
the modulation frequency must be much larger than the 
frequency of tunneling oscillations, u> S> 2tt/T(uj). 

Now we calculate the time evolution using the two-site 
Bose-Hubbard Hamiltonian with J(t) of Eq. ([4]). Here 
we take |JVl,JVr) — |0,JV) as an initial condition. In 
Fig. f2J we show the time evolution of (AN) /N [panel 

(a) ] and its fluctuation a an = y (AN 2 ) — (AN) 2 [panel 

(b) ] for N = 5, U/J = 4, A = 0.1, and Tiuj/J = 15 
}21| as an example. Resulting dynamics has a periodic 
nature with a period T ~ 88To, which is much shorter 
than that without modulation evaluated by Eqs. (j2J and 
©: T = 1228.8T . Here T = 2tt/w = tt^/J, where 
ujq = 2J/H and To are the oscillation frequency and the 
period of the tunneling for U = 0. We also see that 
a an almost reaches a an — N — 5 when (AN)/N = 
corresponding to a NOON state. Figure [^c) shows the 
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FIG. 2: (Color online) Time evolution of (AN)/N (a) and its 
fluctuation ctajv (b) for JV = 5, U/J = 4, A = 0.1, and 
Tiuj/J — 15. The tunneling period is T ~ 88Tq. Panel 
(c) shows the population of each component of |AJV) at 
t = 285.27T . 
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FIG. 3: (Color online) The tunneling period T as a function 
of cj for JV = 5, U/J = 4, and A — 0.1. There is a drastic 
reduction of T in a wide range around huj/J ~ 16 (red arrow). 
The dashed lines show T without modulation. Very narrow 
resonances in the low-o; region of Too/ J ^ 10 are not shown. 
The inset shows for JV = 4, U/J = 5, and A = 0.3. 



snapshot of the population of each component of |AJV) 
at t = 285.27T ; here, we have an almost perfect NOON 
state of cf an = 4.994. 

In Fig. |3l we show the resulting tunneling period T as 
a function of w for the same parameters as in Fig. [2] The 
resonance condition §5§ gives Tluj/J = 16. We note that 
there is a dramatic reduction of T around huj/J = 16 
with a wide width in ui (red arrow) [22l | . Let us first dis- 
cuss the high- and low-w regimes. In the high-w region 
of Too/ J > 30, tunneling period does not depend on uj 
any more and is almost the same as in the case without 



3 




FIG. 4: (Color online) Eigenvalues 4>\an\ of the Floquet oper- 
ator as functions of Tilj/J for the same parameters as in Fig. [3] 
TV = 5, U/J = 4, and A = 0.1. The red points (they look like 
lines except for the low-w region) show 05 and two sets of the 
red lines almost overlap. Resonance occurs at u where 4>5 (red 
lines) crosses the other 4*\/\n\\ i.e., at hu/J ~ 11.3, 14.0, 15.5- 
16,22.6, and 28.0. 
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FIG. 5: (Color online) The tunneling period T around 
hjLoj J = 16 as a function of tko/J for various values of A. The 
other parameters are the same as in Figs. [5] and O TV = 5 
and U/J = 4. From the lower to the higher lines, A — 0.5 
(dotted line), 0.1 (solid line), 0.05 (dashed line), and 0.01 
(dashed-dotted line). 



modulation. This is simply because hu>/ J > 30 is bigger 
than any other energy scales in this problem and thus the 
system is insensible to the modulation of huj/ J > 30. In 
the low u) region of hui/J < 10 (but ui being much larger 
than the frequency of the tunneling oscillation) , u> depen- 
dence of T is weak (except for very narrow resonances, 
which are not shown in this figure). Since 2it/uj -C T 
even though uj is small, the tunneling period can be 
evaluated by the time average of Eq. $2$ with replacing 
J by J(t) of Eq. Q: T lQW u = 2n/ (AE NiJ(t) ) t with 

(AE N>m ) t = (2J(i)TV(J(t)/[/TV) A '- 1 ) 4 /3(TV,TV) = 

2JNn- N+1 f3(N,N){(l + Asmut) N ) t . This replace- 
ment of J N by (J N (t)) t yields an effective increase of 
the energy splitting and thus a reduction of T. 

Next, we discuss the narrow resonances observed in 
Fig. [3] For this purpose, Floquet theory analysis is useful 
(e.g., Refs. [20j,[23j). The Floquet operator F is a map- 
ping between the state at to and the state after one mod- 
ulation period at t + 2it/uj: \^(t Q + 2tt/w)) = F\Hf(t )}. 
Here we get F as follows. Starting from each state of 
the basis set {|ATV);ATV = —TV, —TV + 2, TV}, we fol- 
low the time evolution for one modulation period, 2tt/u!. 
Each resulting state forms a column of the matrix of F. 

In Fig. 01 we show the eigenvalues of F as functions of 
Tllo/ J for the same parameters as in Fig. [3] We denote 
the eigenvalue of F as 4>\an\ whose eigenstate has maxi- 
mum amplitude at ±ATV with even or odd parity in the 
Fock space labeled by ATV (However, there are signifi- 
cant populations in the other components. They become 
larger when two eigenvalues have a crossing in Fig. UJ). 
In this figure, the red points (they look like lines except 
for the Iow-oj region) show <j)§ and two sets of the red lines 
are so close that they appear as a single line. By com- 
paring with Fig. |31 we observe that the resonances occur 



when </> 5 (red lines) crosses the other </>|aat| (black lines): 
the narrow resonances occur at fnu/J ~ 11.3,14.0,22.6, 
and 28.0, where <j)§ crosses <j>i, and the wide one occurs at 
Tiu)/ J ~ 15.5-16, where (f>5 almost simultaneously crosses 
two sets of 03. 

To understand the difference between the wide and 
narrow resonances, we also show the tunneling period 
T around huj/J — 16 for various A in Fig [5] With 
decreasing A, the resonance becomes narrower and fi- 
nally it separates into two narrow resonances (see the 
case of A = 0.01). Note that the separation of these 
two resonances is tiu/ J ~ 0.75, corresponding to the en- 
ergy splitting Ai?3 between the two states of ATV = ±3. 
This fact shows that, to obtain the wide resonance as 
in the case of A = 0.5 and 0.1 (and 0.05) in Fig. [SJ it 
is necessary to couple the two states of |ATV| = TV — 2, 
whose energies arc different by AT5 , at_ 2 , to the states of 
| ATV | = TV (whose energy splitting is negligibly small). 
Thus, the energy scale of the amplitude of the modu- 
lation in the hopping term of Eq. ([1} should be com- 
parable to or larger than the energy splitting of the 
states of | ATV | = TV - 2: 2 J ATV > AE N _ 2 - From 
Eqs. © and ©, we then obtain A > AE N ^ 2 /2JN = 
N-^J/U)"- 3 ^- 1)(TV- 2)/(TV- 3)!. Concerning this 
point, we have some remarks. For TV = 3, this condition 
reads A > 2/3, which might be too large to be realized 
in a real double- well potential. For TV = 2, there is only 
one state of ATV = TV — 2 = 0, and thus we never have 
a wide resonance. On the other hand, for larger TV, fur- 
ther reduction of T is possible by employing additional 
modulations corresponding to wide resonances between 
the states of | ATV| = TV - 2 and TV - 4, TV - 4 and TV - 6, 
etc., whose energy splittings can be very small for larger 
TV [see Eq. ©] 24]. This "scalability" may be useful for 
creating NOON states with larger TV. 

Note also that a tilt V tilt = AVAN/2 of the double- 
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well potential, which may be caused by the imperfec- 
tion of the double well or by the superimposed trap- 
ping, suppresses the tunneling [1, [25|, [26[ . To obtain the 
NOON state using the present scheme, the tilt should 
be smaller than the energy splitting for |AiV| = N — 2: 
AVN « AE N -2- 

It is useful to know the width V of the wide res- 
onance, which enables us to evaluate the formation 
time for a given w. We numerically obtain V by the 
Lorentzian fit of the following form: T — (Tjy — 
T N - 2 ) (1 - T 2 [(u; - cj ) 2 + T 2 ]- 1 ) + Tjv-2. This formula 
is constrained to reproduce T = T^^ 2 at u = wo [I?! 
and T = TV in the limit of u> — > oo. We find that T 
is almost proportional to A and U, respectively. From 
the results of N < 6, we obtain the scaling relation as 
r~ 0.49AU(N- l)(N-2), 

Finally, we discuss the feasibility of creating larger 
NOON states in current experiments 13]. In the present 
situation with N > 2, most serious decoherence process is 
due to three-body losses. Let us estimate the decoherence 
time assuming the Gaussian wave function in each well, 
V>(r) = (N/^l 2 d\d z fl 2 exp [-(r\/2d\) - {z 2 /2d 2 z )] , 
where we take the direction of the double well in the 
z direction and d± and d z are the transverse and the lon- 
gitudinal oscillator lengths, respectively. With the three- 
body rate constant K 3 , the decoherence time T3 is given 
by rg- 1 = K 3 Jd 3 r - (v^tt)- 3 ^^ 3 ^!^)" 1 - 

(At:/3) 3 K 3 N 3 s ± s 1 z /2 {\ 2 l \ 4 z )- 1 . Here X z and A_l are the 
wavelengths of the lasers in the longitudinal and trans- 
verse directions, respectively, and s± and s z are the lon- 
gitudinal and transverse lattice heights, respectively, in 
units of the recoil energy in the longitudinal direction 




In the experiment of Ref. [lj|, X z = 765 nm, A_l = 843 
nm, s z ~ 10 (for U/J = 5L s± = 33, and K 3 = 
5.8 x 1CT 42 Hz m 6 for 87 Rb 29], and we obtain r 3 ~ 16 
ms for N = 4. Using the present scheme, the tunneling 
period for N = 4 and U/J = 5 30] can be reduced to 
< 20T (see inset of Fig. [3]), corresponding to < 14.4 ms 
in this experiment. Thus, the system safely undergoes 
one period of the oscillation in the decoherence time. By 
the same technique used in Ref. (l~3| . one can measure 
(AN). Its oscillatory behavior is a firm evidence for the 
coherency of the system and rules out the possibility of 
the 50:50 mixture of \N, 0} and |0, N). From the oscil- 
lation period of (AN), one can also reject the possibil- 
ity of \N/2, N/2) at (AN) = 0. The combination of 
the above two evidences establishes the realization of a 
NOON state. 

In conclusion, we have found the higher-order co- 
tunneling can be assisted by sinusoidally modulating the 
tunneling matrix element. Using this mechanism, we 
have proposed an efficient scheme for periodically cre- 
ating almost perfect NOON states in two-state Bose sys- 
tems. This scheme is scalable in a sense that, for larger 
number of particles, further reduction is possible using 
more than one frequency. With this scheme, a NOON 
state of four particles can be realized in current experi- 
ments of cold Bose gases in a double-well potential. 
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